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Abstract 

We extend to regular QCD the derivation of a confining qq Bethe-Salpeter 
equation previously given for the simplest model of scalar QCD in which 
quarks are treated as spinless particles. We start from the same assumptions 
on the Wilson loop integral already adopted in the derivation of a semirela- 
tivistic heavy quark potential. We show that, by standard approximations, an 
effective meson squared mass operator can be obtained from our BS kernel and 
that, from this, by ^ expansion the corresponding Wilson loop potential can 
be reobtained, spin-dependent and velocity-dependent terms included. We 
also show that, on the contrary, neglecting spin-dependent terms, relativistic 
fiux tube model is reproduced. 
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I. INTRODUCTION 



In a preceding paper [Q] we have derived a confining Bethe-Salpeter equation for tlie 
simplified model of a scalar QCD. In that we started from the same assumptions already 
used in the derivation of the semirelativistic potential for a heavy quark-antiquark system 
0], 0. The basic object was the Wilson loop integral 

W = hTTFrexpig{<( dxf'AA) . (1.1) 
3 Jr 

where as usual the loop T is supposed made by a quark world line (Ti), an antiquark world 
line (r2) followed in the reverse direction and closed by two straight lines connecting the 
initial and the final points of the two world lines {yi,y2 and xi,X2)', Afj_{x) denotes a colour 
matrix of the form ^4^(2;) = ^A^A"; Pr prescribes the ordering along the loop and Tr denotes 
the trace of the above matrices; the expectation value stands for the functional integration 
on the gauge field alone. 

The basic assumption was 

ilniy = z(lniy)p,rt + o^^min (1.2) 

(Iniy)pert being the perturbative contribution to \iaW and Smm the minimum area enclosed 
by r. Furthermore, we used for S^i^ the straight line approximation, consisting in replacing 
S'min with the surface spanned by the straight lines connecting equal time points on the 
quark and the antiquark worldlines. In practice we wrote 

^min=/ dtr ds[l-{s-^ + {l-s)-^f]-2, (1.3 

Jti Jo at at 



where t stands for the ordinary time in the center of mass frame, zi(t) and Z2(t) for the quark 
and the antiquark positions at the time t, denotes the transverse velocity {6'^'' — ^-^)-^ 
of the particle and r(t) the relative position zi(t) — Z2(t). 

In ref.[l] an essential tool was the Feynman-Schwinger path-integral representation for 
the spinless one-particle propagator in an external field. 
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In this paper we want to extend the derivation to the case of the regular QCD with the 
lagrangian 

L = Y: Mh'D^ - mf)i^f - jF.^F^'' + Lgf (1-4) 
/=i ^ 

(where — d^ — igA^ and Lgf is the gauge fixing term), in which the quarks are fermions 
and have spin and / is the flavour index. 

To do this we find convenient to work in the second order formahsm. 

As usual the gauge invariant quark-antiquark Green function is given by 

Gf (xi,X2,yi,y2) = ^(0|TV'2^(a:2)C/(x2,xi)V'i(xi)V^i(yi)C/(yi,y2)V^2(y2)|0) = 

= ^-Ty{U{x2, yr, A)U{y,, y2)C-'S2{y2, x^, A)C) (1.5) 

where c denotes the charge-conjugate fields, C is the charge-conjugation matrix, U the 
path-ordered gauge string 

U{h, a) = P6„ exp i^g dx^ A^{x)^ (1.6) 

(the integration path being the straight line joining a to 6), 5*1 and 5*2 the quark propagators 
in the external gauge field A^^ and obviously 

/P[A]M,(A)/[A]e-^[^l 

S[Al\ being the pure gauge field action and Mf{A) the determinant resulting from the explicit 
integration on the fermionic fields (which however in practice we set equal to 1 in the adopted 
approximation) . 

The propagators Si and 5*2 are supposed to be defined by the equation 

(Z7^D^ - m)S{x, y; A) = 5\x - y) (1.8) 

and the appropriate boundary conditions. 
Alternatively we can set 

S{x, y; A) = {lYD, + m) A-(x, y; A) (1.9) 



and have 

{D,D^ + m'-^g a^''F^,)A'^{x, y; A) = -5\x - y) (1.10) 
[a'^" — Then, taking into account gauge invariance we can write 

Gf{xi,X2]yuy2) = {nidxin + mi){i-f2dx^i, + ■m2)H^{xi,X2;yi,y2) (l-H) 

with 

H4{xuX2;yuy2) = -^IV([/(x2, xi)A^(xi, yi; l/2)A5(x2, ^2; -i)) (1.12) 

and the tilde denotes transposition on the colour indices alone. What we shall show 
is that the "second order" Green function Hi{xi, X2] yi, ^2) satisfies a Bethe-Salpeter type 
nonhomogeneous equation of the form 

Hi{xi, X2\ yi, ^2) = H2{xi -yi)H2{x2 -y2) - i j d%d%d\d'^r]2 

H2(xi - ^1) H2{x2 - 6) hi^i, 6; Vi, V2) H4{r)i, 772; yi, ^2) (1-13) 

where H2 stands for a kind of colour independent one particle dressed propagator and I4 
denotes the appropriate kernel which is obtained as an expansion in the strong coupling 
constant fts = f;;: and in the string tension a (better in aa^, a being a characteristic length, 
typically the radius of the particular bound state) . 

At the lowest order in ag and aa^ we can write in the momentum space 

Hpi^ P2, Pi, P2) = 4ert(pi, P2; Pi, P2) + -^conf(pi, P2; Pi, P2) (1-14) 

{p'l + P'2 = Pi + P2) with 

(1- 

and 



/conf = J dhe'^ ''J{v, qi, ga) 



:i.l6) 



with 



■Jir, qi, q2] 



2ar 



qio + q20 



+ -; r (arcsm h arcsm ■ 



IqtI 



qio 



q2o 



+ 



^ <y\''q2Qqiur^ , ^ (yl''qwq2y , 

-ZO" = + zcr , = + . . . 



r^qlo + Qt 



:i.l7) 



In (1.15)-(1.17) we have set 
p'l + Pi 

qi = — ^ — ' 



^2 



P2 + P2 



Q=p'i-Pi=P2- p'2 



(1.18) 



Dpa{Q) denotes the gluon free propagator and the center of mass system is understood 
qi = — q2 = q, Q't = {^^'^ ~ f^'-f'^)q^. 

Eqs. (|1.13|) - (|1.17| ) are the basic resuhs of this paper. 

Notice that Eq. ( |1.15 ) corresponds to the usual ladder approximation in this second 
order formalism, while ( |1.17 ) reduces to Eq. (1.8) of |l| when the spin dependent terms 
are neglected. Notice also that instead of ( 1.13|) one could have written the homogenous 
equation 



$p(A;') = -i 



(2vr) 



r^2(r7iP + k')H2{v2P - k')I{k\ k; P)$p(A;) 



(1.19) 
P 



which is more appropriate for the bound state problem. In this equation rjj = ^^^^^ , 
denotes the total momentum pi +p2, k the relative momentum 772^1 — V1P2 {qj = VjP + 
and in the CM frame q = ), ^p{k) is the ordinary Bethe-Salpeter wave function in the 
momentum space. 



From (|1.19| ) by replacing H2{p) with the free propagator 



and performing an ap- 



propriate instantaneous approximation on I [consisting in setting Qq = 0, qjc 



w'.+Wj 



or 



Pjo = PjO 



w'.+Wj 



or fcn 



kn 



= ^2 



uj^wi _ W2+W2 ^ _^ _^ _^ y^2); with 



k^' ] one can obtain an effective square mass operator for the 



mesons (in the CM frame P = 0, P = (m^, 0)) 



= M^ + U (1.20) 



with Mq = W mf + + W m2 + and 



The quadratic form of Eq.( |1.2Cl|) , obviously derives from the second order character of the 
formahsm we have used. It should be mentioned that for light mesons this form seems 
phenomenologically favoured with respect to the linear one. 
In more usual terms one can also write 

M = Mq + V (1.22) 

with 

(k'l^lk) = -^^i_/i„st(k', k) + . . . (1.23) 

l^^J 4:JWiW2WiW2 

where the dots stand for higher order terms in as and aa^ and kinematical factors equal to 1 

2 

on the energy shell have been neglected. In the limit of small ^ the potential V as given by 
( |1.20| )-( p!.22| ) reproduces the semirelativistic potential of ref. 0,0]. Similarly, if we neglect 



in V the spin dependent terms and the coulombian one, we reobtain the hamiltonian of 
the relativistic flux tube model |^ with an appropriate ordering prescription 0. However 
we have not yet completely understood the relation between the spin dependent terms we 
obtain and those appearing in the relativistic flux tube model with "fermionic ends" recently 
proposed [H. 



Finally we want to mention that a result directly in hamiltonian form (|1.22|) but strictly 
related to our one has been obtained by Simonov and collaborators 

The paper is organized in this way. In Sect. II we discuss the Feynman-Schwinger 
representation for the one quark propagator in an external field, in Sect. Ill and IV we 
study the corresponding representation for and derive the BS equation for such quantity, 
in Sect. V we introduce the effective mass operator and consider its semirelativistic limit 
and its relation with the flux tube model. Finally in Sect. VI we summarize the results and 
make some additional remarks. 



II. THE FEYNMAN-SCHWINGER REPRESENTATION 



The solution of Eq.( |l.lCl| ) can be expressed in terms of path integral as (Feynman- 
Schwinger representation) 

i is 1 

A^{x, y; A) = -- dsexp -{-D^D^ - + -ga^^"" F^,) 

j roc rx 1 O 

= --/ ds VzF^yT^yexpt dT{--{m^ + z^)+gA,{z)zP + ^yF^,{z)} (2.1) 

Z Jo Jy Jo Z 4 

where the path integral is understood to be extended over all paths z^ = z^(r) connecting 
y with X and expressed in terms of a parameter r with < r < s, z stands for the 
"functional measure" is assumed to be defined as 

Vz={-^f''d'z,...d'zN-i, (2.2) 
2me 

Pxy and Txy prescribe the ordering along the path from right to left respectively of the colour 
matrices and of the spin matrices. 

On the other side, it is well known that, as a consequence of a variation in the path 
z^(r) z^^ij) + 5z^{t) respecting the extreme points, one finds 



S{F^yexpig rfri;^(r)A^(2;)} 



= ig 6S^''{z{T))P,y{ - F^Mr)) expzg rfr'i^(r')A^(^(r'))} (2.3) 
with 6S'^'^{z) = ^{dz^6z'^ — dz^5z^). So one can naturally write 

T..exp(-1 dra^'js^)) {^^v '9 [ dr' z^{r')AMr'))) 

= T^yV^yexvig f^dT[z^{T)A^{z{T)) + -^a^^ F^,{z{t))] (2.4) 

and Eq. ( p.l|) can be rewritten as 

A''{x,y;A) = -- dr V zP ^yT ,yS^ exp t dT[--{m^ + z^) + tg^'^A^iz)] (2.5) 

Z Jo Jy Jo Z 

with 



= exp - - / dra^"-—— (2.6) 



In (2.6) it is understood that z^{t) has to be put equal to z^{t) after the action of Sq^ . 
Alternatively, it is convenient to write z = z + (, assume that Sq acts on ({t) with 63'^'^ (z) = 
lidzi'SC - dz^dC^) and set eventually C = 0. 

III. THE FEYNMAN SCHWINGER REPRESENTATION 



Replacing (^l5| ) in ( |1.12|) we obtain 



2 poo poo pXl pX2 

Hi{xi,X2]yi,y2)= {-f j dsi ds2 Vzi P^2T^iyiT^2?/2'5o''5o' 

Z Jo Jo Jyi Jy2 

exp(— ){ / dTi{ml + zl)+ dT2{ml + zl)} 
I Jo Jo 

i(TrPr expitg){ j( dz^A.iz)}) (3.1) 

where now z = zj = zj + (j on Fi and ^ = 2; on the end lines X1X2 and ?/22/i and the final 
limit (j — > is again understood. 

Then, let us try to be more explicit concerning Eq. (1.2) and (1.3). For the first term in 
(1.2) we have, at the lowest order of perturbation theory, 

i{\nW)pert = iln(-TrPexpii5( i dzf'A^{z))p^,t = -g / dri / dT2D^^{zi - Z2)z^zl^ 
6 Jr 6 Jo Jo 

2 n /"^l /"^l 2 n pS2 

/ dn / dT[D^,{z^ - z[)z'^z[' - -g^ / dT2 / dT^D^,iz2 - z'^)z^z'^' + ... (3.2) 

3 JO JO 3 Jo Jo 

On the other side, for the second one in general we have to write 

x'^ = u'^{t, s) being the equation of the minimal surface with contour F. Let us assume that 
for fixed t and for s varying from to 1, u^{s^t) describes a line connecting a point on the 
quark world line Fi with one on the antiquark world line F2, 

^'^(l,^) = ^r(ri(t)), u^{Q,t) = z^{T2{t)) (3.4) 

Obviously ( ^.3| ) is invariant under reparametrization, so a priori t and s could be everything. 
In particular if Fi and F2 never go backwards in time, t can be choosen as the ordinary time 
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u'^[s,t) = t. For the moment let us assume that this is the case. Then Ti{t) and T2{t) are 
specified by the equation 

z'M) = zl{T,) (3.5) 

and we can set 



Obviously L cannot depend only on on the extremal points -Zi(ti) and -22(^2) but has to 
depend even on the shape of the world lines at least in a neighbourhood of such points. So, we 
can think of it as a function of all derivatives in ti and T2 and write L = L{zi, Z2, ii, Z2, ■ ■ ■)■ 
Finally (3.3) can be rewritten as 

5'min = J dz^ J dz^S{z'^ - Z^)L{zi, Z2, Zi,Z2,...)=J J dT2S{z^ - Z^)z^Z^L{Zi, Z2, Zi, Z2, . . 

(3.7) 

and in principle this expression can be considered a good approximation even if the world 
lines contain pieces going backwards in time. In fact, in such a case if we fix e.g. ri, ( p.5| ) has 
more than one solution in T2 and if Fi and F2 are not too much irregular in space ( otherwise 
S'min is large and the weight of the loop is small) the minimal surface can be reconstructed 
as the algebraic sum of various pieces of surface. 

In the straight line approximation we must choose 

u^{s,t) = t 

u\s,t) = sz\{n{t)) + (1 - s)z\{T2{t)) (3.8) 

and we have 

z^zlL = a\z[ - c?s{ifo4o - ("SZiT^o + (1 - ■5)z2T^io)^}' (3-9) 

Jo 

which introduced in (3.7) becomes equivalent to Eq.(1.3). The important point concerning 
(3.7) with (3.9) is that it has the same general form as (3.2). However we stress that the 



approximation (3.8) must be performed only after thai the application of the operators Sq^ 
and Sq^ has been performed. 

Substituting (|3l^) and (3.7) in (p.l[) we obtain 



roo roo rxi nX2 

H4{xi,X2;yi,y2) = (t:)^ dsi ds2 Vzi Vz2TxiyiTx2y2<So^<So^ 

Z Jo JO Jyi Jy2 

X fsi 1 rS2 

expi{--/ dTi{ml + zl)-- dT2{ml + z^) + 
Z Jo Z Jo 



/ dn [ dT[Df,^{zi - z[)z'^zl' +'^g^ [ dT2 [ dT2D^i,{z2 - z'^jz^z. 
o Jo Jo 6 Jo JO 



Si rs2 



^1 \ i^' 
2 

dTi I dT2E{zi,Z2,Zi,Z2,...)}, (3.10) 

Jo 
where we have set 

4 

E{zi, Z2, zi, Z2---) = -g'^D^y{zi - Z2)Z^Z2 + 

+(y5{ziQ - Z2o)zioZ2oL{zi, zs, h . . .). (3.11) 

Now let us denote the quantity in curly bracket in (3.10) by S4 and perform a Legendre 
transformation by introducing the momenta pj^ = ( in this the various quantities Zj, 

i 

z j ^ z j ^ • • • are assumed to be treated as independent) 

• ^^af'^'n (- -<\-'^ I A , dE{zuz'2.Zi,Z2...) 
= z^i + -g dT^D^^[zi - z^)z^ + dT2 ^rfi 

• ^^2r',/n f '\->^ I r A , dE{zi,z'^.zi,Z2...) 

Pf,2 = z^2 + dT^D^^[z2 - Z2)Z2 + dr^ . (3.12) 

Eq.( p.l2D cannot be inverted in a closed form with respect to ii and Z2. However, we can 
do this by an expansion in = |^ and aa^ and at the lowest order we have 

4„2r'^_/n .'^^'u ^^,dE{zuz'2,Pi,p2-..) ^ 



^1 = Pi - / dT[D^,y{zi - z[)p'^ - / c/r2 
3 JO Jo 

•M M 4 2 , /n (- -i\ IV n ^ ,dE{zi,z'2,p[,p2...) ^ 

Z2=P2-i^9 / dT2D^^{z2 - z^jp^ - dr^ — + ... (3.13) 

3 JO JO 0P2 

with 

P^=P^ + C^- (3.14) 

In conclusion we find 
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- X2y2 



]^ roo poo rxi PX2 

i^4(Xi,X2,?/l,y2) = (tt) / dSi dS2 VZiVpi T>Z2'Dp2T^^y^T^ 

I Jo Jo Jyi Jy2 

dTiKi + dT2K2- dri dr2E(^i, ^2, Pi, P2, •■•) + •■■ (3-15) 



where 



1 2 /"^j 

= -Pj ■ Zj + 2^p'j- ^]) + dTjD^^{zj - z'j)p'^p 



(3.16) 



includes the self-interacting term. Notice that now in Sq^ it must be understood 5S^^{zj] 



Eq. ( p.l5|) is the starting point for the derivation of our Bethe-Salpeter equation. 



IV. THE HOMOGENEOUS BETHE SALPETER EQUATION 



In Eq. (|3.15| ) we proceed along the same line followed in Ref. [0]. 
Applying to the interaction term E the identity 



exp / dTA{T) = l+ / rfrA(r)exp( / rfr'v4(r')) 
Jo JO Jo 



(4.1) 



we have 



poo poo pXl pX2 

H^{xi,X2;yi,y2) = (-)^ / dsi / ds2 / VziVpi / Vz2'Dp2 

Z JO Jo Jyi Jy2 

( rsi rS2 rsi r-S2 

TxiyiT^2y2^o'^0^\^^P'^[j^ '^^^^^"^i) ^^2i^2]-«y^ dT2E{Zi,Z2,Pl,P2- ■ ■) 

Usi rS2 rri rS2 ~\ 

dnKi + J^ dT2K2-J^ dT[ dT!^E{z[,z'^,p[,p'^,...)^ (4.2) 



Now, using the method of Ref. |]l| and having in mind ( p.ll| ), ( p.6|) and (|3.1CI| ) one finds 
(see Appendix for details) 
S 



6S^"'{z 



rSi rS2 

— dT[ dT^E{z[, 4, p'l, P2, . . .) 
ijJo Jo 



S2 



dr2 -g'^{dyD^„{zi - 4) - d^D^^{zi - z'^)pl + 



V\v[ziy^ — zL) — pi^i{zii, — Z2^) 



10 -L3 

+a5{ziQ - Z20, , ^ 

V(p?o-pD(zi-Z2)2 + (pi-(zi-z'2))2 

and a similar result, with a minus sign of difference in front, for the derivative in ^ 



(4.3) 
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Furthermore 

5St"'{zi)SSP''{z[) Jo Jo SSi"'{z2)6Sp''{z'2) Jo Jo 

but 

A2 PSl 1-32 4 

(4.5) 



5Sp^''^{zi)5Sp^''^{z2) 



Then, taking into account the relation 

n=0 

and specifically ( [1.4|) , we have 



i?e-^ = E-T[A[A...[Ai?]...]], (4.6) 



5o^^5o^^ dn r dr^Eiz,, Z2,Pi,P2, . . = 
Jo Jo 

j dTi j dT2E{Zi,Z2,pl,p2,---) = RiZi,Z2,Pl,P2) (4.7) 



with 



R — -Rpert + -Rconf (4.8) 



4 , 

^pert = --g'^[Dp^{zi- Z2)p1p2 



3' 

1 



--crn^pi, - S^MDpA^i - ^2)PI - -,<{^ld2u - 6:d2,)D,^{z, - Z2)p{ 



+^<yr'or\KA^. - dtA,,){6,.d2.2 - d:A,,)D,„{z^ - Z2)] (4.9) 



16 



and 



-Rconf= 0-5(Zio - ^20)||Zl " Z2I ds\JplQplQ - [spiTP20 + (1 " s)P2TPio]^ 

-4^20^1 . , / , + -^P10CJ'2 1 / 2 /• (^-^^^ 

^ |Zi-Z2|-y/p|o-pfT |Zl -Z2IVP2O -P2T ^ 

Finally setting 
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H2{x-y) = ^ ds VzVpS'oexpi drK (4.11) 

2 Jo Jy Jo 

Eq.(1.18) can be written 

H4{xi,X2; Vi, 2/2) = ^^2(2:1 - yi)H2{x2 - ^2) + 

fOO fOO fXi nX2 PSl rS'z 



^ /-OO roo rxi nX2 rsi rS2 

-- / dsi / ds2 / VziVpi / p2;2pp2T^is/iT^'2y2 / / dT2R{zi, Z2,pi,P2) 
4 Jo Jo J1/1 Jj/2 JO Jo 

S'o'S^'expz{ £dT[K[ + fj dr'2K'2-i £ dT[ rfr^E(z;, 4,pi,p2 • • •)}• (4-12) 



At this point it is necessary to take explicitely into account the discrete form of ( [4.121) . 
If we take 

P exp [ig j( dz^A^iz)] = P [] U{z^, z^_,) = P exp ig ^(z^ - ^^i)A^( ^" + ^"'^ ) (4.13) 
(as required by a gauge invariant definition of the integral on the gluon field) we have 



00 00 

6 



Hi{xi, X2\ 2/1, 2/2) = ii2{xx - y\)B.2{x2 - 2/2) - 

^ iV;^=0 Ar2=0 

(27r)^i+^2 y d'^Vwd'^^W ■ ■ ■ d'^^lN-ld'^PlNi J d^P2ld'^Z2i ■ . ■ d"^ Z2N2-ld^P2N2-l^ xiyi"^ X2y2 

EV^ T)/ZlRl + ^IR-L-l Z2R2+Z2R2-I ^ 
^( 7. > 7y ,PlR,P2R) 

Rl=l iJ2 = l ^ ^ 

2 iVj 

'^o'^o"^ exp i{ Z Z [ - - 2;jn-i) + e{p% - m|) - 

j=ln=l 

-gf ^ L ^ Y ^^On'VjA 

n'=l 

Rl-1 N2 ^ , ^ r -U r 

2 j-,/^ln "T ^ln-1 ^2n "T ^2n.-l n1 . n 

Z ^( 2 ' 2 'Pi"'P2n,...)| (4.14) 

ni=l 712=1 

If we neglect in the exponent the "non planar" terms YIULr^+i Y,n'=i ^iJ.'^Pjn^P'jn'. 
„2^=/?2+i ^f^'^PinP2n2' corresponding in the continuous to the quantity 

r dr'^ r dr^D^y (4 - (4. 15) 

Jr, Jo 



and 



dT[ r dT^D,yiz[ - (4.16) 

Jo Jr2 
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Eg. ( [4. 141) can be written 



• oo oo oo 

i 



H4{xuX2,yuy2) = H2{xi-yi)H2{x2-y2) - js'^ J2 J2 

^ iJi=l R2=l A^i=Ri+l A'2=R2+l 

— ^ / d'^ ZiR^d^piR^<f ZiR^-ij—^ I d'^Z2R2d^p2R2d^ZR2-i / VziVpi / Vz2Vp: 



(27r)^ J ^i^i"- /^i-ni"' ^2 



'2 

Vz2Vp2 

yi •'y2 

2 Ni ^ Ni 

T^l^i«T^2^2«'S^l542 exp J2 [-PjniZjn - Zjn-l) + £ip% - "^i) - of H ^M^PinPin' 

j=ln=i?j+l n'=i?j+l 

exp [ - tl^PjnizjR^ - ZjR^_,)]R{ , ,Pir,P2r) ■ 

j 

2 Rj-1 Rj-1 

T^z^nyiT^Z2Ry2<So'So^ expz{ ^ ^ [ - Pjn{Zj„ - Zjn-l) + £{P% - m^j) - Y D^^PjnPjn'] + 

j=l n=l n'=l 

+ Y: Y ^( '^"^/^"-\ Pln„ J>2„2)} (4.17) 

ni=l 712=1 

which, going back to the continuous corresponds to Eq.( |1.13D with 



/(^i,6,?/i,^2) = -4z j ^2iTY^ ^i ^ 2 ^^ '^^'^^-'^^P^^"^-'^'^^^ -^i)^i + (6-^2)/^2]} 

(4.18) 



In conclusion, taking the Fourier transform 

{2ny6{p[+p',-p^p2)i{p[,p'2■,Pl,P2) = -4z / d%d% J d\d% 

we obtain (1.14)- ([n7|) . 



V. EFFECTIVE MASS OPERATOR 



As we mentioned, for bound states Eq. (1.13) can be replaced by 

$p(A;') = -tJ -^H,{r]iP + k')H2iv2P - k')i{k\ k; P)<I>p(fc) (5.1) 



(2vr) 

(P = P1+P2, k = r]2Pi — 'riip2, k' = T]2p'i — 'riiP2 and in the center of mass frame P = 0, Pq 
V^, Pi = -P2 = k, p'l = -P2 = k'. 
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Let us recall the definitionof the instantaneous kernel given in Sec.l and consider the 
approximation consisting in replacing in (5.1) I{k',k;P) by /inst(k, k') and in substituting 
H2j{p) with the free propagator ^2^^- Further let us introduce the reduced wave function 



(^p(k') 



2wi(k)«;2(k') 



/'OO 
dk'o^p{k'). (5.2) 
-OO 



and integrate over ko and kg using 

1 



{kg + r/ims)2 — k^' — ml + ie {—kg + ri2mBy — k^' — m2 + ie 
w'^W2 rriQ — [w'l + W2y 



we have 
(«;i(k')+w;2(k'))Vm,(k') + 



+ / (2^^ 2^,(k0^2(kO -^'-^(^^'^^)^ 2^,(k)^2(k) V^"-B(k)-m.^^,(k) (5.4) 

from which Eq.(1.21) immediately follows. 

The linear potential of Eq.(1.22) is then given by 

(k'|V|k) = {k'\U\k) + ... (5.5) 

^ ^ ^ ' W[ + W'2 + Wi+W2^ ^ ^ ' ^ ^ 

from which (neglecting kinematical factors becoming equal to one for |k| = |k'|) we obtain 
Eq.(1.23). By performing a ^ expansion on Eq. (1.23) we find the qq potential at ^ order 



3 '27r2Q2 7r2Q4 3 27r2mim2Q2' ' 

-o^"s k — - ^ k + £hki^^{<yi + (^2 k ^ k 

3 2mi r"* 2m2 r"* 3 2mim2 2 r"* 

+ |3^ ^ k) (7^(72 + -CTi ■ CT2 

3 mim2 r'^ 3mim2 3 

-^(^ + ^-^)(k1q^r|k) 

mf m2 17111712 



■77£/»fei^T^ — — + — )(k — k) + i^(k + 

2 2 mf m2 r 2 mi r m2 r 



with q = Q = k' - k. 
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Passing to the coordinate representation we may also write 

.2 )^'^'}w 



V= -— + ar + {— 5'^" + 

6 r 6 mim2 2r r 



lyt /X ^ 

hk , ' ' \„h„k^ 



3 2m^ 



1 ai ■ r 1 ^2 ■ r 4 

+ -7^ (o-i + 0-2) • (r X q) 



1 7 - 2m2 ' 3 2mim2 
1 as 3(cri ■ r)((T2 • r) 0-1-0-2, 4 27r 



6mim2 r° 
a ^ 1 
--(— + — 

^2 77111712 



3 mim2 3 



>i ■(T2)5^(r) 



){qTr}w 



o- 0-1 



CTZ 1 «! ■ r 1 0:2 ■ r 



2 mf m2 r 2 mi r m2 r 



(5.7) 



where now q stands for the momentum operator. Now, by performing a Foldy-Wouthuysen 
tranformation with generator 



S 



-ai ■ q - - — a2 ■ q 



2mi 2m2 

we end up with the ^ potential which coincides with the Wilson loop potential 
4 a 



(5.8) 



V 



3 r 



+ or 



im2 13 r J w 



2mim2 

2 1 1 
- ^^{^'"PjTlw - 7 Wr^pi^ ■ P2t}w 



1 



4as 
3 r 



+ ar + 



-^Si ■ (r X pi) - ^82 ■ (r X P2) 



1 Mas a 

2 V3 r'^ r 
1 Aa 

+ -^[82 ■ (r X pi) - Si ■ (r X P2)] + 



m 



rrin 



+ 



17111712 3 



+ 



1 



mim2 3 



-(Si-r)(S2-r)-Si-S2 



+ ^53(^)8^.82} , (5.9) 



with r = (r/r) and the symbol { }w stands for the Weyl ordering prescription among 
momentum and position variables. 



VI. CONCLUSIONS 

In conclusion, under the assumption (1.2) and (1.3) for the evaluation of the Wilson 
loop integral, we have derived a quark-antiquark Bethe-Salpeter (BS) equation from QCD, 
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extending a preceding result obtained for spinless quarks. The assumptions are the same 
previously used for the derivation of a semirelativistic heavy quark potential and the tech- 
nique is strictly similar. The kernel is constructed as an expansion in ag and aa^ and at the 
lowest order is given by equations (1.14)- (1.18). 

The BS equation that has been obtained is a second order one, analogous in some way to 
the iterated Dirac equation. Correspondently, by instantaneous approximation, an effective 
square mass operator can be obtained from (1.19) which is given by (1.20) and (1.21). 

At the lowest order in and aa^ even a linear mass operator can be written with a 
potential V given by (1.23). Neglecting the spin dependent terms in V the hamiltonian 
for the relativistic flux tube model comes out. On the contrary by a ^ expansion and 
an appropriate Foldy-Wouthuysen transformation the ordinary semirelativistic potential is 
reobtained. 

In equation (1.13) or (1.19) a colour independent dressed quark propagator appears which 
is defined by equations (4.11) and (3.16). Notice that only the perturbative expansion gives 
contribution to this quantity. 

Few additional remarks are in order. 

First of all, notice that the result does not depend strictly from equation (1.3) or (3.9) 
but from the possibility of writing the interaction term as an integral on the world lines of 
the quark and the antiquark, as evidenced in (3.10). Multiple integrations of the same type 
would be admissible, as it occurres for the perturbative contribution, but dependence of the 
integrand on higher derivatives in the parameters ti and r2 would not enable to carry on 
the argument. We have no actual justification that i In is in general of the desired form, 
we observe however that this quantity is obviously independent of the parametrization. For 
an example of inclusion of higher order perturbative terms see Ref . [1] . 

A second point concerns the significance of the lowest order BS kernel we have derived. 
As the analysis in terms of potentials show, the inclusion of terms in ctg is essential for an 
understanding of the fine and the hypcrfinc structure. For what concerns the importance of 
cr^ contributions an indication can be obtained considering the corresponding terms in the 
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relativistic tube flux model. Neglecting the coulombic terms and in the equal mass case the 
cm. hamiltonian for such model at the order can be written 



ar V \/ + . | qrp | / + q^ 



-arcsm 



qxl i/m^ + q2 y + q' 



+ 



+ 



^2^2 ^2 _|_ q2 |- ^ ffi^ -(- q2 _ |q^| /m^ + q^ 



- arcsm - 



2 



(6.1) 



16qT -\/m2 + q2 L |qrp| A/m^ + q^ ym^ + q^- 

To better appreciate the relative magnitude of the two potential terms let us consider e.g. 
the case of small qx (small angular momentum) in which the above equation becomes simply 

2 2 

i?cm = 2v^m2 + q2 + ar + — (6.2) 
* 16v^ + q 

Then, taking into account that a ~ l/(crm)^, g ~ 1/a, and assuming typically a = 

0.17GeV^, mu = 0.35GeV, vric = 1.7GeV, mb = 5GeV we find that the last term in 

(6.2) is of the order of the 5%, 0.8%, 0.2% of the preceding one for the uii, cc, bb systems 

respectively. This would correspond to contributions to the mass of the meson of about 20, 

2, 0.2 MeV. The inclusion of the coulombic term would reduce a and improve the result. In 

the uu case e.g. it would amount to cut the above contribution by a factor 2. Therefore 

only in this last case the terms would be of any significance. 

Finally let us come to the problem of the type of confinement, which has been largely 

discussed in the literature. By this terminology it is usually meant the tentative assumption 

of a BS (first order) confining kernel of the instantaneous form 

Jeonf = (27r)3r4ij, (6.3) 

or even the covariant counterpart of it 

/conf = -(27r)3r-^ij, (6.4) 

where F is a combination of Dirac matrices. Typically the cases F = 1 (scalar confinement), 
F = 7i72 (vectorial confinement) or a combination of them have been considered. 

Eq. ( |6.4D is immediately ruled out by the fact that, even if formally it corresponds to 
( |6.3|) (by instantaneous approximation), actually, due to the strong infrared singularity, it 



gives results very different from ( |6.3|) . As well known, Eq. ( |6.3| ) with F = 1 was motivated 
by the fact that it reproduces the static potential ar and the spin dependent potential as 
obtained in the Wilson loop context. This choice, however, gets both into phenomenological 
and theoretical difficulties: 

1. it gives a firts order velocity dependent relativistic correction to the potential 
which differs from the Wilson loop one and does not seem to agree with the 
heavy meson data 0, 

2. it does not reproduce straight line Regge trajectories ||10|,3- 



Complementary objections can be moved to ( |6.3| ) with F = 75^72- 

On the contrary, even if we have not yet attempted calculations directly with the kernel 
established in this paper, very encouraging results have been obtained in the context of 
the relativistic flux tube model of the dual QCD and of the effective relativistic 
hamiltonian |^, formalisms that are all strictly related to our one. Therefore the complicated 
momentum dependence appearing in (1.16)-(1.17) seems essential to understand both the 
light and the heavy meson phenomenology. 

APPENDIX A: APPENDIX 

We want to prove Eq.(4.3). 

Let us first consider the confinement part and rewrite Eq. (3.3) as 

^min = r dt dsS{u) (Al) 
Jt; Jo 



with 



S(u) 



~dt~dt] \ ds~ds] ^ \ dt ds ] 



(A2) 



Being x'^ = u'^{s,t), the equation of the minimal surface enclosed by the loop must be 
the solution of the Euler equations 
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d dS 



d dS 



(A3) 



satisfying the contour conditions u^{l,t) = 2;f (ri(t)), ■u^(0,t) = Z2{T2{t)). Then, considering 
an infinitesimal variation of the world line of the quark 1, Zi{t) — > Zi{t) + Szi{t), even 
u^{s, t) must change, u^^{s, t) — u^{s, t) + 6u^^{s, t) and one has 



SSrr 



Jti Jo 



dt 



ds 
dS 



dS d . ,, dS d ^ 



[dm 



5u^ 



dt 



(A4) 



where Szi{t) is assumed to vanish out of a small neighbourhood of a specific value of t. 
Finally taking into account that 



one obtains 



rt( 1 

SSuiin — dt -T—: 

= I /*' dt {dz'^5z'[ - dz'[Sz'^) 
2 Jti 



ds 



1 



ds 



Sz'i 



X 



ds 



' du 



Zlu 



' duy 
ds 



/ du' 



and more explicitely 



Then, in the straight line approximation we have 

dUa 



ds 



— Zin — Z2fj, — 



and 



dS„ 



dS^^-{zl) [-zlr^j^{zi-rY]^ 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 
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Using Eg. ([A9|) (having substituted the velocities with the momenta), Eqs.(3.11) and (3.6), 
we obtain the second term in (4.3). 

Let us come to the perturbative part. Consider a variation zi zi + 6zi, then 

Si J dn J dT2z{Df„{zi - Z2)z1 = 

= J dn J dT2[5z{Dp„{zi - Z2) + z{5zidyDp„{zi - Z2)]z2 = 

= J dS"" J dT2[d,Dp„{zi - Z2) - dpDUzi - Z2)\zl (AlO) 



and so 



J dn J dr2p{DpM - Z2)pl = J dT2i6;di, - 5Pdip)Dp,{zi - Z2)pl (All) 



5^"'{zl 

and then we recover the first term in (4.3) 
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